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摘要 GWAC 的 超大 视 场 每 天 可 观测 到 数 百 个 流星 体 , 识别 这 些 流星 体 可 为 流星 的 科学 研究 提供 重要 的 
科学 数据 。 文章 针对 GWAC 这 一 类 超大 视 场 测 光 巡 天 系统 在 流星 识别 时 ， 遇 到 的 不 能 有 效 区 分 流星 与 
其 他 移动 目标 的 挑战 ， 设 计 和 实现 了 一 个 流星 候选 体 识别 算法 ， 该 算法 主要 包含 流星 轨迹 识别 和 光 交 
曲线 形态 分 析 两 个 部 分 。 识 别 算法 对 Mini-GWAC 约 两 个 月 的 图 像 进行 处 理 ， 提 取 10.9 万 个 移动 目标 
ik, HP 90% 以 上 属于 非 流星 目标 。 分 析 午 夜 时 间 段 内 流星 目标 高 斯 拟 合 曲线 的 aR, AMAR 
分 单 峰 流星 目标 的 光 变 曲线 波峰 随 图 像 像素 位 置 呈 慢 速 变化 趋势 。 综 合流 星 的 单 帧 特性 、 光 变 曲 线 的 
单 峰 结构 特征 和 光 变 曲线 的 慢 速 变化 特性 进行 过 滤 ， 最 终 得 到 4.1% 的 高 精度 流星 候选 体 。 经 过 人 工 检 
查 确认 ，4.1% 的 流星 候选 体 中 有 85%~87.3% 的 目标 符合 流星 的 形态 及 亮度 特征 。 

关键 词 流星 ; 拉 长 轨迹 识别 ， 光 变 曲 线 ; 形态 分 析 ; GWAC 
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1. 引言 

地 基 广 角 相 机 阵 (Ground Wide Angle Cameras， 简 称 GWAC， 详 见 第 2 节 ) 是 中 法 合作 天 文 卫星 SVOM 
的 地 基 相 关 观 测 设备 ， 主 要 科学 目标 是 探测 伽 玛 射线 暴 、 引 力 波 事 件 等 爆发 时 对 应 的 电磁 辐射 对 应 体 。 除 了 探测 
{FSS Be aR Ey, GWAC 约 5000 平方 度 的 超大 视 场 同时 能 够 探测 到 大 量 的 移动 目标 ， 如 流星 、 近 地 目标 、 
小 行星 、 彗 星 、 飞 机 、 衬 间 和 碎片、 卫星 等 。 这 些 移动 目标 会 影响 伽 玛 暴 等 暂 现 源 的 探测 效率 ， 识 别 并 过 滤 这 些 移 
= 动 目标 对 GWAC 系统 具有 重要 意义 。 同 时 ， 国 际 上 的 许多 组 织 在 进行 各 类 移动 目标 的 研究 ， 如 国际 天 文联 盟 的 
一 ”流星 数据 中 心 致力 于 流星 数据 的 收集 和 流星 特性 的 研究 中 , GWAC 系统 将 会 贡献 出 大 量 的 流星 数据 。 本 文 主要 针 
对 GWAC 中 的 流星 目标 设计 识别 算法 。 
际 上 有 许多 流星 研究 项 目 : NASA 资助 建立 的 Cameras for Allsky Meteor Surveillance project (CAMS) 项 
BI 包含 3 个 站 点 共 60 个 录像 相机 ， 其 目标 是 系统 研究 流星 雨 的 特性 及 来 源 ; 日 本 的 SonotaCo" 4; 25 个 站 点 
共 100 个 录像 相机 , 已 经 向 国际 天 文联 盟 的 流星 数据 中 心 贡献 了 大 量 的 流星 两 数据; 欧洲 的 European viDeo Meteor 
Observation Network (Edmond) 流星 数据 库 中 汇聚 了 欧洲 多 个 流星 观测 网 的 数据 ， 克 罗 地 亚 的 Croatian Meteor 
Network (CMN) 基于 廉价 的 监控 相机 和 普通 台式 机 构建 面向 科普 教育 的 流星 观测 系统 ， 加 拿 大 Western 大 学 
的 Southern Ontario Meteor Network (SOMN) "主要 关注 厘米 级 流星 体 的 探测 ， 俄 罗斯 的 天 文科 学 望远镜 阵 
Mini-MegaTORTORA (MMT) 与 GWAC 系统 有 着 类 似 的 科学 目标 ， 即 探测 短 时 标 光 学 暂 现 源 ,MMT 包含 900 
平方 度 浅 视 场 和 100 平方 度 的 深 视 场 ， 白 光 0.1 秒 曝光 极限 星 等 约 11V，1.5 年 累计 观测 到 约 9 HREL 

现 有 项 目 使 用 的 流星 识别 软件 如 MeteorScan!'®!, MetRec!'!!, UFO Capture!’?!, AIM-IT"?!, Gural ^il vida"?! 
等 ， 主 要 面向 流星 巡天 项 目 进行 设计 ， 其 图 像 的 主要 特点 是 : D 极限 星 等 低 ( 约 SV 星 等 )， 视 场 中 背景 恒星 的 
数目 少 ， 识 别 流星 时 干扰 源 少 ; 2) 帧 频 快 〈 每 秒 曝 光 几 十 到 上 百 帧 )， 拍 摄 的 流星 体 轨 迹 出 现在 多 帧 图 像 中 ， 每 
帧 图 像 中 包含 一 段 轨 迹 ， 通 过 多 段 轨迹 可 计算 出 流星 体 的 角速度 : 3) 多 站 点 观测 ， 流 星 巡 天 项 目 通常 包含 少 则 
两 个 多 则 几 十 个 观测 站 点 ， 联 合 两 个 及 以 上 站 点 的 观测 数据 ， 可 以 计算 出 流星 体 的 位 置 、 高 度 和 线 速度 。 流 星 巡 
天 项 目 中 流星 识别 的 过 程 为 : 1) 提取 多 帧 图 像 中 的 轨迹 线段 ， 关 联 多 个 轨迹 线段 为 一 个 目标 ; 20 通过 多 站 点 数 
据 ， 计 算 目 标的 位 置 ， 3) 通过 多 个 轨迹 线段 的 位 置 计算 目标 的 速度 : 4) 通过 目标 的 高 度 和 速度 信息 对 流星 候选 
体 进行 过 滤 。 王 恩 旺 59 使 用 改进 的 帧 差 法 对 图 像 中 拉 长 轨迹 很 短 的 空间 运动 目标 进行 检测 。 

GWAC 图 像 的 特点 是 : 1) 极限 星 等 高 (GWAC 为 16V 星 等 ，Mini-GWAC 为 12.5V 星 等 )， 图 像 中 会 有 大 量 
的 背景 恒星 及 移动 目标 ; 20 "UAE CEE 10 秒 ， 读 出 5 秒 )， 拍 摄 的 流星 出 现在 单 帧 图 像 中 ， 通 常 表现 为 一 段 
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完整 的 拉 长 轨迹 〈 见 2.2 节 )， 无 法 计算 流星 体 的 速度 : 3) 目前 GWAC 项 目 仅 有 一 个 观测 站 点 ， 无 法 计算 流星 
体 的 位 置 和 高 度 。 由 于 GWAC 的 图 像 和 图 像 中 流星 的 轨迹 形态 与 现 有 流星 巡天 项 目的 有 着 本 质 的 区 别 ， 因 此 现 
有 流星 识别 软件 或 算法 不 能 直接 应 用 于 GWAC 的 流星 探测 。 

际 上 与 GWAC 有 着 类 似 科 学 目标 的 大 视 场 测 光 巡天 项 目 有 MMT" RIT Pi of the Sky ， 这 两 个 项 目 中 也 没 
有 可 借鉴 的 经 验 : 1) MMT 的 曝光 时 间 为 0.1 秒 ， 其 流星 识别 算法 与 流星 巡天 项 目的 类 似 ; 2) Pi of the Sky 没有 
发 表 流 星体 识别 的 相关 研究 成 果 。 


- 


^ 


图 1 Mini-GWAC 图 像 中 各 种 类 别 的 目标 ， 其 中 上 面部 分 是 原始 图 像 ， 下 面部 分 是 相 邻 图 像 
相 减 后 的 结果 。(1) 为 流星 候选 体 ; (2~6) 为 不 同类 别 的 移动 目标 ; (7) 为 望远镜 视 场 边缘 被 遮挡 后 
产生 的 假 目 标 ; (8) 为 流星 的 尾 焰 ; (9) 为 热 像 列 
Fig.1 All kind of objects in Mini-GWAC, top part of image is original image, bottom part of image is residual image of adjacent image. 
(1) is meteor candidate, (2~6) are movingobjects, (7) is fake object, (8) is wake flame of meteor, (9) is hot column 


图 REREH, MAI TREAT, SLEPA LAME 7 
Fig.2 The center part of meteor is wider and brighter than both ends 

GWAC 系统 的 超大 视 场 每 天 不 仅 能 观测 到 数 百 个 流星 体 ， 同 时 也 能 观测 到 几 千 个 其 他 类 型 的 移动 目标 ， 如 
流星 、 飞 机 、 卫 星 等 。 如 图 1 所 示 ， 为 Mini-GWAC 〈 见 2.1 节 ) 图 像 中 各 种 不 同类 别 的 移动 目标 ， 除 (1) 是 流星 
候选 体 之 外 ， 其 他 为 不 同类 别 的 干扰 源 。 与 其 他 移动 目标 相 比 ， 流 星 轨 迹 的 典型 特征 为 1》 持续 时 间 短 ， 普 人 
单 帧 出 现 ， 2) 其 轨迹 在 形态 上 中 间 宽 两 端 罕 ， 在 亮度 上 中 间 亮 两 端 暗 〈 火 流星 例外 )， 如 图 2 所 示 。 

GWAC 系统 的 流星 体 识别 过 程 面 临 的 挑战 是 : 在 无 法 依靠 速度 和 高 度 信息 准确 区 分 流星 体 和 其 他 移动 目标 
的 情况 下 ， 如 何 从 大 量 的 移动 目标 中 找到 流星 体 。 文 章 将 针对 GWAC 的 系统 特性 及 其 图 像 中 流星 的 轨迹 特征 设 
计 流 星 识别 算法 : 1) 针对 流星 的 拉 长 轨迹 和 单 帧 特征 设计 流星 轨迹 识别 算法 ; 2) 针对 流星 轨迹 的 亮度 和 形态 特 
征 ， 对 流星 候选 体 的 形态 特征 进行 分 析 ， 以 进一步 提高 流星 的 准确 率 。 文 章 在 第 二 部 分 简要 介绍 GWAC 系统 和 
分 析 GWAC 系统 中 移动 目标 的 轨迹 特性 ， 第 三 部 分 对 流星 轨迹 识别 算法 和 形态 分 析 算 法 进行 详细 描述 ， 第 四 部 
分 简要 介绍 了 算法 的 实现 细节 , 第 五 部 分 对 部 分 移动 目标 轨迹 的 形态 进行 统计 分 析 , 最 后 一 部 分 总 结 了 算法 的 优 
缺点 ， 及 对 下 一 步 工作 的 展望 。 
2. GWAC 系统 简介 
2.1. GWAC 观测 系统 


图 3 GWAC 系统 ， 左 边 为 Mini-GWAC, i27; GWAC 
Fig.3 Image on the left is Mini-GWAC, and on the right is GWAC 


GWAC 系统 的 建设 包含 两 个 阶段 , 分别 是 一 期 工程 Mini-GWAC 实验 系统 和 二 期 工程 GWAC, 如 图 3 所 示 。 

1) Mini-GWAC 的 组 成 : 共 包 含 12 台 6 厘米 的 宽 视 场 望远镜 ,每 个 望远镜 配备 一 台 3KX3K 的 CCD 相机 ， 
综合 视 场 约 5000 平方 度 ， 恒 星 探测 极限 星 等 约 12.5V (10 秒 钟 、56o、 无 月 夜 );， 于 2014 年 10 HÆR, 
并 于 2015 年 10 月 正式 开始 观测 。 

2) GWAC 的 组 成 : 共 包 含 40 台 有 效 口径 为 18 厘米 的 宽 视 场 望远镜 ， 每 个 望远镜 配备 一 台 4kX4k 高 性 能 
CCD 相机 ， 综 合 视 场 约 5000 平方 度 ， 恒 星 探测 极限 星 等 约 16.0V C10 秒 钟 、5o、 无 月 夜 )。 现 在 已 完 
成 部 分 望远镜 的 调试 ， 并 在 进行 试 观测 。 

Mini-GWAC 已 经 积累 了 大 量 的 观测 数据 ， 文 章 使 用 Mini-GWAC 的 数据 进行 流星 识别 算法 的 研究 ， 并 在 后 
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续 应 用 于 GWAC. 
2.2. Mini-GWAC 系统 移动 目标 轨迹 特性 分 析 


Mini-GWAC 图 像 中 的 移动 目标 主要 包含 3 类 : 飞机 、 流 星 和 卫星 《空间 碎片 ) ， 下 面 将 简要 分 析 这 三 类 目 
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标的 特性 。 
A 飞机、 流星 和 卫星 的 高 度 和 线 速度 
Table 1 The height and speed of airplane, meteor and satellite 
最 低 高 度 (km) ”最 高 高 度 (km) ”最 小 线 速 度 (km/s) KERERE (km/s) 
飞机 6 12.6 0.08 0.3 
流星 80 120 11 72 
卫星 200 36000 3 7.78 
1) ”高 度 和 速度 特性 


如 表 1 所 示 ， 飞 机 、 流 星 和 卫星 三 类 移动 目标 的 高 度 和 速度 区 间 有 着 明显 的 差异 ,计算 出 高 度 或 速度 中 的 一 
即 可 区 分 移动 目标 的 类 别 。GWAC 系统 目前 为 单 站 点 观测 ， 无 法 计算 移动 目标 的 高 度 或 速度 ， 因 此 无 
法 利用 高 度 或 速度 特性 对 移动 目标 进行 区 分 。 


项 参数 ， 


2) 


时 间 特 性 


MMT 与 Mini-GWAC 使 用 相同 型 号 的 镜头 ， 


0.1 秒 到 2.5 秒 外 。 照 此 估算 ， 流 星 在 Mini-GWAC 


其 在 1.5 年 中 累计 观测 到 约 9 万 颗 流星 ， 这 些 流星 持续 时 间 为 


的 图 像 〈10 秒 曝光 5 秒 读 出 〉 中 只 能 单 帧 出 现 。 飞 机 、 卫 星 


等 移动 目标 在 飞 出 图 像 边 缘 前 会 持续 存在 ， 因 此 可 以 在 多 帧 中 连续 出 现 。 


3) 


轨迹 长 度 特性 


图 4 所 示 为 对 表 1 中 三 类 典型 移动 目标 在 Mini-GWAC Eti 10 秒 曝 光 图 像 中 轨迹 长 度 的 估算 。 可 见 在 不 同 
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4) 


plane height(km) 
6 8 10 


速度 和 高 度 时 , 这 三 类 移动 目标 在 图 像 中 的 轨迹 长 度 存 在 交集 ， 因此 无 法 直接 通过 轨迹 长 度 对 这 三 类 移动 目标 进 
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图 4 Mini-GWAC 的 一 帧 10 秒 曝光 图 像 中 ， 流 星 、 飞 机 和 卫星 的 轨迹 长 度 估算 ， 
其 中 流星 的 长 度 按 0.5 秒 持续 时 间 估算 ， 飞 机 和 卫星 的 长 度 按 10 秒 持续 时 间 估算 
左 图 为 流星 与 飞机 的 轨迹 长 度 对 比 ， 右 图 为 流星 和 卫星 的 轨迹 长 度 对 比 
Fig.4 The estimation of track length of airplane, meteor and satellite in a 10 seconds exposure image of Mini-GWAC, 
the left picture shows the track length comparison of the meteor and the plane, and the right picture shows the 
track length comparison of the meteor and the satellite 


轨迹 形态 特性 


Subasinghel 1 对 CAMS 双 站 点 观测 的 891 颗 流 星 进 行 了 光 变 曲线 和 形状 分 析 : 1) 光 变 曲线 分 析 : 光 变 曲线 
的 形状 有 67% 的 单 峰 目标 ， 包 括 14% 的 早 峰 类 型 ， 
C110 帧 / 秒 ) 中 流星 为 类 华星 状 。 流 星 在 Mini-GWAC 图 像 中 单 帧 出 现 , 其 轨迹 也 包含 光 变 曲线 和 形状 信息 : 1) 


42% 的 对 称 类 型 ，11% 的 晚 峰 类 型 ，2) 形状 分 析 : 单 帧 图 像 
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受 噪声 的 影响 ， 因 此 本 文 不 对 流星 的 形状 进行 分 析 。 
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曲线 : 从 Mini-GWAC 的 图 像 中 能 够 得 到 流星 随 位 置 变化 的 光 变 曲线 ， 后 文 将 会 进行 详细 介绍 ; 20 形状 特 
fE: Mini-GWAC 图 像 中 流星 的 轨迹 中 间 粗 两 端 逐 渐变 细 ， 其 宽度 仅 为 几 个 像素 到 十 几 个 像素 ， 采 样 不 足 ， 容 易 


3. 流星 候选 体 识别 算法 
3.1. 流星 轨迹 识别 
算法 对 时 空 连续 的 多 张 FITS 图 像 1) 进行 处 理 ， 整 体 流程 如 图 5 所 示 ， 下 面 将 着 重 对 几 个 关键 的 步骤 进行 
介绍 。 
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图 5 流星 轨迹 识别 算法 的 流程 图 
Fig.5 The flow chart of meteor track recognition method 
D 图 像 预 处 理 
连续 的 两 帧 图 像 拍摄 时 间 相 隔 15 秒 ， 图 像 之 间 seeing 等 环境 因素 的 变化 非常 微弱 ， 可 直接 对 连续 两 幅 图 像 
进行 相 减 操作 : 


变 亮 ， 移 动 目 


停止 循环 过 程 。 


S(k)- I(k)- I(k-1) 为 图 像 顺序 编号 ， 取 值 为 2.…N) (1) 
对 相 邻 的 两 幅 图 像 相 减 , 避免 了 天 文 数 据 处 理 流程 中 对 暗 场 和 本 底 的 操作 , 同时 也 不 用 对 图 像 进行 去 背景 的 
操作 。 原 始 图 像 相 邻 素 点 之 间 是 有 关联 的 ， 且 相 邻 图 像 之 间 有 细微 的 局 部 噪声 起 伏 〈 如 热 噪 声 、 大 气 扰动 等 ) ， 
为 了 降低 环境 因素 对 关键 信息 的 影响 ， 可 对 相 减 后 的 图 像 SEIT AAE. SOIA EREA Rk), HA 
理 可 表示 为 : 


相 减 后 图 像 SHERUAR 3 部 分 内 


A: 背景 〈 像 素 值 接 近 零 ) ， 残 差 信 号 〈 局 部 信号 起 伏 ) ; 真实 信号 〈 星 
标 等 ) 。 图 像 二 值 化 过 程 包括 : 


W(k)= R(k)"? S(k) (2) 


1) 对 W(k)f H1 30 RUS Sete, MABE E Albee D; 


2) 初始 化 阔 值 设 为 T=E+2D， 使 用 阀 值 卫 对 S(K) 进 行 二 值 化 操作 得 到 二 值 化 图 像 B(k): 


255, W(k) 2T 


BUE i , Wk) <T 


(3) 


3) 如 果 公式 (5) 中 霍 夫 变换 操作 没有 得 到 线段 集 ， 且 T>100, WIES) AS BEL T: 


循环 每 次 将 T 的 值 减 5， 然 后 对 图 


2) ”移动 目标 轨迹 识别 


T=E+2D-N*5(N 为 正 整 数 ) (4) 
像 进 行 二 值 化 和 线段 提取 操作 。 如 果 T«100 或 者 公式 (5) 找 到 线段 集 ， 则 


文章 选用 概率 霍 夫 变换 (progressive probabilistic hough transform, PPHT) 03 来 对 二 值 化 图 像 进行 线段 探测 。 


相 比 霍 夫 变换 ， 概 率 霍 夫 变换 计算 速度 快 ， 且 能 直接 算出 线段 的 两 个 端点 坐标 。 概 率 霍 夫 变 换 后 ， 得 到 候选 线段 


集 : 


L(k)-PPHT(B(k)) (5) 

流星 在 Mini-GWAC 图 像 中 的 轨迹 宽度 约 几 个 像素 到 十 几 个 像素 。 因 此 霍 夫 变换 后 ， 一 个 流星 轨迹 有 可 能 被 
识别 为 多 条 线段 。 为 去 除 元 余 线段 ， 文 章 对 候选 线段 集 L(K) 进 行 融 合 ,将 相 邻 的 多 条 线段 合并 为 一 条 线段 。 文 章 
综合 采用 线段 中 心 距离 、 线 段 斜 率 和 图 像 中 心 到 线段 的 距离 三 个 参数 来 判断 线段 的 相似 性 ,对 这 三 个 参数 小 于 指 
定数 值 的 线段 合并 为 一 条 线段 Lo Lo 代表 一 个 移动 目标 候选 体 Oo， 对 线段 集 L(k) 合 并 后 得 到 移动 目标 候选 集 
O(k). 

3) ”多 帧 过 滤 

Mini-GWAC 图 像 中 不 仅 包含 流星 ， 还 包含 其 他 的 移动 目标 〈 如 飞机 、 卫 星 等 ) 。 流 星 与 其 他 移动 目标 的 区 
别 是 : 流星 仅 在 单 帧 图 像 中 出 现 〈 流 星 的 持续 时 间 为 0.1~2.5 秒 ，Mini-GWAC 单 帧 图 像 曝光 时 间 10 秒 ， 读 出 时 
间 5 秒 ) ， 而 其 他 移动 目标 可 能 会 在 多 帧 图 像 中 出 现 。 因 此 ， 对 合并 后 的 移动 目标 候选 集 O( 刀 进行 前 后 帧 关联 ， 
可 以 过 滤 一 部 分 非 流 星 目 标 。 
3.2. 流星 光 变 曲线 形态 分 析 

移动 目标 光 变 曲线 从 形态 上 可 分 为 无 峰 目 标 、 单 峰 目 标 和 多 峰 目 标 三 类 ， 如 图 6 rz. Subasinghe" 分 析 的 
流星 样本 中 有 67 多 的 流星 是 单 峰 目标 ， 且 其 主要 针对 单 峰 目 标 进 行 光 变 曲线 的 形态 分 析 。 因 此 文章 尝试 对 单 峰 
目标 的 形态 进行 分 析 ， 以 期 进一步 过 滤 非 流星 目标 。 由 于 流星 在 Mini-GWAC 图 像 中 主要 为 单 帧 出 现 ， 无 法 得 到 
流星 随时 间 变 化 的 光 变 曲线 ， 但 是 可 以 通过 流星 在 图 像 上 的 拉 长 轨迹 得 到 流星 随 图 像 像素 位 置 变化 的 光 变 曲线 ， 
后 文 简称 光 变 曲线 。 下 面 将 简要 介绍 光 变 曲线 的 提取 及 单 峰 目 标 形态 分 析 的 关键 步骤 。 

LS <r — À 


图 6 三 类 移动 目标 ， 左 图 为 无 峰 目 标 ， 中 图 为 单 峰 目标 ， 右 图 为 多 峰 目标 


Fig.6 Three kinds of moving objects, left is no-peak object, center is single-peak object, right is multi-peak object 
1) ” 光 变 曲线 的 提取 
2x Bt Lo 5j i EVER X 轴 的 夹 角 为 8， 对 Sk) ET HEH 0,， 使 目标 长 轴 在 图 像 中 处 于 水 平 位 置 ， 对 目标 进 
行 矩形 区 域 裁 图 得 到 如 图 6 所 示 的 窗口 图 Subo。 对 Subo 沿 图 像 Y 轴 方 向 对 X 轴 方向 的 流量 强度 进行 求 和 , 得 到 
目标 亮度 随 x 的 变化 曲线 : 


Flux, = Xy Subo(x, y) (6) 

2) 峰值 计算 及 多 峰 目 标 过 滤 

峰值 计算 过 程 参考 SExtractor29 中 的 多 颗 相 邻 恒星 目标 的 解析 过 程 。 对 光 变 曲线 Flux, W X 轴 平 行 的 直线 
y=Yo, Yo 从 光 变 曲线 Fluxo 的 最 大 值 开 始 ， 依 次 以 Step, 递减 ， 直 到 最 小 值 。 每 次 递减 后 计算 直线 与 光 变 曲线 的 
交叉 部 分 , 交叉 部 分 的 个 数 即 为 峰值 个 数 peakNum 。 如 果 峰 值 个 数 peakNum 等 于 0, 则 代表 该 目标 为 无 峰 目 标 ; 
如 果 峰 值 个 数 peakNum 等 于 1， 则 代表 该 目标 为 单 峰 目 标 ， 如 果 峰 值 个 数 peakNum 大 于 1， 则 代表 该 目标 为 多 
峰 目标 。 

3) ” 光 变 曲线 形态 分 析 


文章 采用 高 斯 函数 拟 合 的 方式 对 光 变 曲线 进行 形状 的 量化 , 多 个 高 斯 函数 车 加 能 更 准确 的 描述 目标 的 光 变 曲 
线 ， 这 里 至 多 使 用 两 个 高 斯 函数 对 光 变 曲线 进行 拟 合 。 光 变 曲 线 拟 合 函 数 如 下 : 
Flux(k) = Y. Ai * ii ,N-[12] (7) 


j 拟 合 后 高 斯 函数 的 参数 (4 、xoz 和 分 别 代表 第 i 个 拟 合 曲线 的 高 度 、 中 值 和 峰值 锐 度 ) 作为 对 目标 光 变 
曲线 形态 的 量化 参数 。 文 章 将 在 第 5 节 对 目标 的 光 变 曲线 形态 进行 详细 分 析 与 讨论 。 
4， 程 序 实 现 

文章 使 用 Python 和 OpenCV 对 算法 进行 实现 ， 在 i 4770k CPU 上 测试 ， 每 幅 图 像 的 处 理 时 间 不 超过 1 秒 ， 
满足 GWAC 系统 的 实时 探测 的 需求 。 算 法 实现 的 关键 信息 如 下 : 

1) ”异常 图 像 过 滤 : Mini-GWAC 的 单 图 像 视 场 覆盖 天 区 范围 大 ， 易 受 环境 因素 的 影响 。 文 章 通过 相 令 图像 
相 减 后 残 差 的 方差 和 有 效 像素 百分比 这 两 个 参数 对 图 像 进 行 过滤 。 
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2) S AER AXE: 移动 目标 轨迹 经 过 二 值 化 之 后 ， 有 可 能 会 形成 多 个 线段 。 通 过 分 析 Mini-GWAC K 
像 中 流星 的 特征 ， 选 取 概 率 霍 夫 变 换 测 参数 为 p=1 像素 ，0=1 度 ， 组 成 轨迹 的 子 线段 最 小 长 度 为 10 像素 ， 子 线 
段 间 最 大 间隔 为 5 像素 ， 轨 迹 的 最 短 长 度 为 50 像素 。 

3) ”线段 融合 参数 选择 : 对 图 像 中 多 条 线段 ， 如 果 中 心 距离 小 于 50， 且 线段 斜率 差 小 于 $， 且 线段 与 图 像 中 
心 的 距离 差 小 于 150， 则 这 些 线段 合并 为 一 条 线段 。 

4) “ 相 邻 帧 线段 集 栈 构造 : 对 连续 多 帧 图 像 产 生 的 线段 集 L(k) 绥 存在 一 个 堆栈 中 , 用 来 过 滤 多 帧 出 现 的 移动 
目标 。 对 相 邻 的 两 个 线段 集 ， 对 一 个 线段 集中 的 所 有 线段 都 与 另 一 个 线段 集中 的 线段 进行 比较 ， 如 果 两 个 线段 的 
斜率 之 差 绝 对 值 小 于 5， 则 认为 这 两 个 线段 是 同一 个 目标 的 线段 ， 即 该 目标 是 多 帧 出 现 的 移动 目标 。 

5) ”峰值 数量 peakNum 的 计算 : 首先 对 流量 曲线 进行 中 值 滤波 , 滤波 器 的 kernelSize 选取 目标 长 度 的 10%。 
峰值 计算 过 程 Stepy=Fluxmaw10。 为 过 滤 噪 声带 来 的 毛刺 ， 限 定 相 邻 两 个 峰值 的 最 大 值 之 差 应 大 于 Stepy 

6) ” 光 变 曲线 归 一 化 : 将 光 变 曲线 的 长 度 和 高 度 (亮度 ) 都 归 一 化 到 [0,1] 之 间 。 对 于 所 有 光 变 曲线 ， 记 最 大 
KEX Lnae， 记 最 大 高 度 为 Horx; 将 所 有 光 变 曲线 等 比 放 大 到 长 度 为 Las: 对 所 有 光 变 曲线 的 长 度 除 以 Lmax 
完成 长 度 归 一 化 ， 对 所 有 交 变 曲线 的 高 度 除 以 Hiae， 完 成 高 度 的 归 一 化 。 

7) ” 拟 合 参数 融合 : 对 于 单 峰 目 标的 拟 合 ， 有 些 目标 O7 使 用 一 个 高 斯 函数 能 有 效 拟 合 ， 有 些 目标 0; 需 要 两 
个 高 斯 函数 才能 有 效 拟 合 。 对 于 01:， 算 法 直接 使 用 拟 合 的 三 个 参数 4、xo 和 a ) 描述 其 形态 。 对 于 Or MEK 
两 个 高 斯 函数 的 参数 进行 融合 ， 以 更 准确 地 描述 流星 曲线 的 形态 : 对 两 个 4 取 最 大 值 ， 对 两 个 xo 取 均值 ， 对 两 
个 x 取 最 小 值 。 

5， 光 变 曲 线形 态 特 征 分 析 
5.1. 数据 统计 
文章 对 Mini-GWAC 的 131 万 幅 图 像 ( 约 两 个 月 的 数据 量 ) 进行 分 析 处 理 ， 结 果 中 各 类 目标 的 统计 如 表 2 所 


D 样本 选择 : 这 里 选择 单 像素 均值 ) 信 噪 比 大 于 等 于 56 的 移动 目标 ， 总 数 约 10.9 万 个 。 

2) EAS WMS: 包含 1.2 万 个 单 帧 出 现 的 目标 和 9.7 万 个 多 帧 出 现 的 目标 。 

3) ” 按 光 变 曲线 峰值 个 数 统计 : 包含 约 2.5 万 个 单 峰 目标 、1.3 万 个 多 峰 目 标 和 7.1 万 个 无 峰 目 标 。 

4) ”流星 候选 体 : 按照 Mini-GWAC 系统 中 流星 的 单 帧 出 现 特性 对 移动 目标 进行 统计 ， 流 星 候选 体 约 占 总 数 
的 10.7%。 如 果 通 过 光 变 曲线 的 单 峰 特性 进一步 过 滤 ， 则 流星 候选 体 的 数据 量 可 降低 到 5.8%。 


表 2 Mini-GWAC 不 同类 别 的 移动 目标 数量 统计 
Table 2 Statistic of different categories of moving objects in Mini-GWAC 


统计 量 数量 (个 ) ”所 占 百分比 
移动 目标 总 数 108764 100% 
NE os ms s a I ade 
多 帧 出 现 数量 97151 89.3% 
—M ng  — Á. —— oC 
"uU 0 0 ÁO M a, 
无 峰 目 标 数 量 70976 65.2% 
Sees e cec ee 
ARES AT O2 的 目标 数量 5588S 
单 帧 出 现 单 峰 峰 值 参数 ar 大 于 02 4453 4.196 
且 拟 合 度 R2 大 于 0.7 的 目标 数量 
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为 量化 光 变 曲 线 的 拟 合 精度 ， 文 章 采用 R-Square(R2) 进 行 评估 。 表 2 中 的 6275 个 单 帧 出 现 单 峰 目 标 光 变 曲 
线 高 斯 拟 合 结果 R 的 直方 图 如 图 7 所 示 。 其 中 有 41.2% 的 目标 R2 大 于 0.9; 有 81.8% 的 目标 R2 KF 07. HF 


数据 处 到 


时 存在 噪声 的 干扰 ， 所 以 选择 较 低 的 R2 un Ce 
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图 7 单 帧 出 现 单 峰 目 标 (6275 个 ) 光 变 曲线 高 斯 拟 合 度 R2 直方 图 
Fig.7 The histogram of R2 of 6275 single-peak objects appeared in one frame 


5.2. 光 变 曲线 形态 分 析 

高 斯 拟 合 参 数 w 代 表 移 动 目标 光 变 曲线 的 波峰 宽度 : w 值 越 小 ， 波 峰 越 窜 ，w 值 越 大 ， 波 峰 越 宽 。 文 章 将 基 
于 KK 对 单 帧 出 现 单 峰 目 标的 光 变 曲线 进行 统计 分 析 。 为 了 找到 能 在 光 变 曲线 形态 上 区 分 不 同 移动 目标 的 参数 ， 
文章 选择 3 个 样本 集 ， 并 对 这 3 个 样本 集 的 波峰 宽度 值 a 进行 直方 图 统计 分 析 ， 这 3 个 样本 集 分 别 是 : 


1) 


2) 


3) 


单 帧 出 现 单 峰 目 标 样本 集 ， 简 称 样本 集 1: 其 中 大 部 分 为 流星 ， 且 包含 少量 的 多 帧 未 匹配 成 功 的 空间 碎 
片 、 卫 星 、 飞 机 等 长 时 间 持 续 存 在 的 移动 目标 。 其 波峰 宽度 值 g 分 布 如 图 8。 


傍晚 或 黎明 期 间 的 单 帧 出 现 单 峰 目 标 样本 集 ， 简 称 样 本 集 22 从 1) 中 选择 晚上 9 点 之 前 和 凌晨 4 点 之 
后 出 现 的 目标 ， 其 中 的 目标 类 别 与 1) 基本 相同 。 其 波峰 宽度 值 a 分 布 如 图 9 所 示 。 


午夜 期 间 的 单 帧 出 现 单 峰 目 标 样本 集 ， 简 称 样本 集 3: 从 1) 中 选择 晚上 10 点 到 凌晨 2 点 之 间 出 现 的 目 
标 ， 该 时 间 段 的 天 空 处 于 地 球 的 阴影 区 ， 除 了 上 自发 光 的 目标 “如 流星 、 飞 机 等 ) ， 很 难 观 测 到 其 他 依靠 
反射 太阳 光 才 能 看 到 的 目标 (如 空间 碎片 、 卫 星 等 ) ; 飞机 等 自发 光 的 目标 受 限 于 速度 , 多 为 多 帧 出 现 ， 
单 帧 出 现 的 概率 非常 低 ， 因 此 该 样本 集中 的 目标 主要 为 流星 ， 其 波峰 宽度 值 gq 分 布 如 图 10. 
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图 8 单 帧 出 现 单 峰 目 标 (6275 个 ) a 参数 的 直方 图 
Fig.8 The histogram of a parameter of 6275 single-peak objects appeared in one frame 
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图 9 晚上 9 点 之 前 和 凌晨 4 点 之 后 ， 单 帧 出 现 单 峰 目 标 (2962 个 ) o 参数 的 直方 图 
Fig.9 Before 9 p.m. and after 4 a.m,the histogram of a parameter of 2962 single-peak objects appeared in one frame 
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图 10 晚上 10 点 到 凌晨 2 点 之 间 ， 单 帧 出 现 单 峰 目 标 (1306 个 ) a 参数 的 直方 图 
Fig.10 Between 10 p.m. and 2 a.m, the histogram of a parameter of 1306 single-peak objects appeared in one frame 


对 这 3 个 样本 集 a 参数 的 直方 图 分 布 曲线 分 析 如 下 : 


1) DIES: 样本 集 1 的 分 布 曲线 有 两 个 峰值 。 这 两 个 峰值 代表 两 类 目标 : 一 类 是 亮度 随 图 像 像素 位 置 快 
速 变化 的 目标 ， 简 称 亮度 快速 变化 目标 ， 其 a 参数 小 于 0.2， 占 比 10%; 一 类 是 亮度 随 图 像 像素 位 置 慢 
速 变化 的 目标 ， 简 称 亮度 慢 速 变化 目标 ， 其 a 参数 大 于 0.2， 占 比 90%. 


2) “分 析 图 9: 样本 集 2 的 分 布 曲线 在 整体 上 与 样本 集 1 基本 一 致 ， 但 样本 集 2 中 的 快速 变化 目标 的 比例 比 
样本 集 1 高 。 
3) ”分 析 图 10: 与 样本 集 1 和 样本 集 2 相 比 ， 样 本 集 3 的 分 布 曲线 在 两 端的 比例 偏 少 ， 其 在 左边 的 快速 变 
化 部 分 也 有 一 个 很 小 的 峰值 , 表明 其 中 快速 变化 目标 所 占 的 比例 非常 少 。 样本 集 3 中 的 移动 目标 主要 为 
流星 ， 因 此 流星 中 快速 变化 的 个 体 非常 少 。 
综合 分 析 图 8、 图 9 和 图 10, 发 现 大 部 分 单 峰值 流星 为 亮度 慢 速 变化 的 目标 , 依 此 条 件 对 样本 集 1 进行 过 滤 ， 
可 以 得 到 5583 的 流星 候选 体 , 约 占 所 有 移动 目标 的 5.1%. 同时 采用 R 大 于 0.7 的 条 件 过 滤 后 得 到 4453 (4.1%) 
个 目标 ， 如 表 2 所 示 。 人 工 对 这 4453 个 目标 进行 检查 ,确认 其 中 有 3785 个 目标 符合 流星 的 形态 及 亮度 特征 ， 占 
比 85%; 有 104 个 目标 处 于 图 像 边缘 ， 图 像 中 仅 包 含 目 标的 部 分 轨迹 ， 这 些 部 分 轨迹 也 符合 流星 的 形态 及 亮度 
特征 ， 占 比 2.3%; FR 564 个 为 低 信 噪 比 目 标 、 天 气 原 因 导 致 的 假 目 标 、 与 流星 形态 类 似 的 天 体 目标 等 ， 占 比 
12.7%。 因 此 算法 的 综合 识别 精度 为 85%~87.3%。 
6. 总 结 与 展望 
GWAC 的 超大 视 场 中 每 天 会 探测 到 大 量 的 流星 ， 有 效 的 识别 出 这 些 流星 能 为 世界 流星 研究 提供 大 量 的 有 效 
数据 。 目 前 在 GWAC 这 一 类 超大 视 场 测 光 巡天 系统 中 进行 流星 识别 ， 遇 到 的 最 大 挑战 是 没有 充足 的 信息 来 有 效 
区 分 流星 与 其 他 移动 目标 。 为 了 尽 可 能 准确 的 识别 出 流星 ， 文 章 设 计 和 实现 了 一 个 面向 GWAC 这 类 系统 的 流星 
轨迹 识别 算法 ， 同 时 对 得 到 的 流星 候选 体 进行 了 光 变 曲线 形态 分 析 ， 结 果 总 结 如 下 : 
1) ”识别 算法 对 Mini-GWAC 约 两 个 月 的 历史 图 像 进行 处 理 ， 从 中 往 选 出 单 像素 信 品 比 大 于 5o 的 移动 目标 
约 10.9 万 个 ; 
2) 通过 流星 在 Mini-GWAC 图 像 中 的 单 帧 出 现 特性 ， 可 过 滤 89.3 多 的 移动 目标 ， 剩 余 10.7% 的 候选 目标 ; 
3) “大 部 分 流星 的 光 变 曲线 为 单 峰 结 构 〈 见 3.2 小 节 ) ， 该 特性 可 进一步 过 滤 候 选 目标 ， 剩 余 5.8% 的 流星 
候选 体 ; 
4) “通过 分 析 3 个 样本 集 光 变 曲线 的 w 参 数 直 方 图 ， 发 现 大 部 分 流星 为 亮度 慢 速 变化 的 目标 。 使 用 该 特性 
和 光 变 曲线 拟 合 度 R2 进一步 对 候选 目标 进行 筛选 ， 最 终 得 到 4.1% 的 流星 候选 体 。 
5) 人 工 对 4.1% 的 流星 候选 体 进行 检查 ， 确 认 其 中 有 85%~87.3% 的 目标 符合 流星 的 形态 及 亮度 特征 。 
文章 中 筛选 流星 使 用 了 光 变 曲线 单 峰 结构 和 亮度 慢 速 变化 这 两 个 特性 , 这 两 个 特征 仅 能 包含 大 部 分 流星 ， 而 
漏 掉 一 些 流星 (如 光 变 曲线 有 多 个 峰值 的 流星 ) 。 但 是 通过 这 两 个 特征 的 过 滤 ， 最 终 得 到 的 4.1% 流 星 候 选 体 具 
有 更 高 的 准确 度 。 
通过 分 析 单 峰 目标 的 实际 图 像 ， 我 们 发 现 一 些 移动 目标 (卫星 、 飞 机 等 ) 的 形状 和 光 变 曲线 形态 与 流星 非常 
相似 ， 如 果 这 些 目标 出 现在 图 像 边 缘 ， 有 很 大 概率 仅 单 帧 出 现 ， 本 文 算 法 很 难过 滤 这 类 非 流星 目标 。 论 文 下 一 步 
将 会 围绕 三 个 方面 对 算法 进行 改进 : 10 增加 对 移动 目标 的 形态 分 析 方 法 ， 如 对 移动 目标 的 形状 进行 分 析 ， 找 到 
多 个 可 量化 的 形态 描述 参数 , 以 期 进一步 提升 流星 的 识别 数量 和 识别 精度 ; 22GW AC 正在 计划 建设 第 二 个 站 点 ， 
两 个 站 点 协同 观测 ， 可 计算 移动 目标 的 高 度 信息 ， 通 过 高 度 信息 可 过 滤 飞 机 和 卫星 等 移动 目标 ; 3) GWAC 正在 


Ur 


ll iioii CMOS 相机 ， 图 像 曝光 时 间 可 达到 0.1 秒 ， 将 极 大 的 增加 移动 目标 的 时 间 分 辨 率 ， 可 计算 移动 目标 


的 速度 信息 。 算 法 中 加 入 高 度 和 速度 信息 后 ， 能 更 准确 的 区 分 流星 与 其 他 移动 目标 。 
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Abstract: With its large field of view, GWAC can record hundreds of meteors every day. These meteors are 
valuable treasures for some meteor research groups. It is therefore very important to accurately find all of these meteors. 
To address the challenge of precisely distinguishing meteors from other elongated objects in a GWAC-like sky survey 
system, we design and implement a meteor candidate recognition algorithm, including the recognizing and morphology 
analysis of the light curves of the meteor candidates. After processing the images of Mini-GWAC taken in two months, 
we detect 109,000 elongated objects in which more than 90 percent of objects are not meteor. By analyzing a parameter 
of gaussian fitting curve of meteor target in midnight, we found that the wave crest of most single-peak meteor target 
changes slowly with the position of image pixel.Among the elongated objects, about 4.1% objects are identified as 
meteors with high confidence, after the filters based upon an existence in a single frame, a single peak in the light curves, 
and a slow variation of the light curves.After manual examination, 85%~87.3% of the 4.1% candidates were consistent 


with the shape and brightness of the meteor. 
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